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Remarks on the calculation of characteristic temperatures 


of cubic crystals from the elastic constants 


By E. TENERz 


Summary 


The Debye characteristic temperature @p for cubic crystals is calculated from the elastic 
constant, using the method first introduced by Bhatia and Tauber (1954). These authors ex- 
pand the mean velocity of the three elastic waves in Kubic Harmonics K,, of which they take 
only three into account. Extending that method by taking five K,,, of which the fifth is deduced 
here, a relative simple expression for 6p is still obtained. It is accurate even for large elastic 
anistrophy. A numerical comparison with Quimby and Sutton’s (1953) method shows a very 
good agreement in the cases where this method is convenient to use. The numerical work is 
almost the same for both the methods in these cases, otherwise much less when Kubic Har- 
monics are used. Values of 0p for a number of substances and temperatures are calculated. The 
new calculation of 9p for KCl at 0° K exhibits good agreement with measurements in the 
T?-range. The fifth K, can be of interest in other calculations. 


Introduction 


In terms of the phase velocities Debye’s characteristic temperature 4p is 
given by the well-known formula: 


h (9S\* 
>= 5 (x) doits 0) 
1 ; 
where I= | Sq p 20 (2) 


Here h and k are Planck’s and Boltzmann’s constants respectively, s denotes 
the number of atoms in a cell, A is the volume of the cell and v;(9¢) are the 
velocities of the three elastic waves for a direction (0¢). The integration in (2) 
is extended over all directions of propagation. 

For a cubic crystal the three velocities, associated with a given propagation 
vector are, according to the theory of elasticity, found from the roots of the 
cubic equation: 

ue —u? + (1— K*) Tu—(1—3 K? +2 K*) 7~=0 (3) 
where 
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Ci, are the elastic constants, 9 is the mass density and 2, y, z are the direction 


cosines of the propagation vector referred to the cubic axis. 


Calculation of (> by expansion of I into Kubic Harmonics 


The Kubic Harmonics K, (x, y,z) were first introduced by von der Lage and 
Bethe [5]. These functions are linear combinations of spherical harmonics so 
constructed as to be invariant under the symmetry operations characteristic of 
a cubic lattice. The problem of constructing them is, firstly, to find combina- 
tions of terms of the form x? y%z'(l=p+q+t) which fulfil the symmetry rela- 
tions and, secondly, to combine them so that they are orthogonal to each com- 
bination with lower 1. 

In von der Lage and Bethe’s paper only the first four K,, are given.. For 
the purpose of this note K, has also been calculated. The expression for the 
first five Kubic Harmonics used in this calculation are: 


K,=1; K,=h,(2'+y'+2-?r); K,=k (ey? 2+3[K)°-17); 
K,=ky (x? +y° +2°—B[K,] 7? — 28 [Ki] — 37°); 


Ky = hey (a +? +2 — 8 [Kg] 9? — 8 [Ka] 1° — 88 [Ky 1 — 89") 


where 7? =2?+y?+z" and km are the normalization constants. The functions in 
square brackets are not normalized, i.e. they represent, Ky, K,/k,, K,/k, and 
Raikes. 

Expanding the integrand in (2) in Kubic Harmonics we get: 


I= | > am KmndQ=42 0%. (4) 
m=0 


An approximate value of «) can now be obtained by using a finite number of 
terms in the expansion, determining the phase velocities for a number of direc- 
tions equal to the number of terms in the expansion and finally solving for the 
coefficient om. This method, first used by Bhatia and Tauber [1], leads to a 
relatively simple expression for J, if only the first few terms are taken into 
account. Bhatia and Tauber, howewer, used only the first three K,,, which 
seems sufficient for crystals where the factor of anisotropy [7 =2 ¢44/(¢,, —¢y9)] 
is not too far from one. In order to calculate 6) for crystals with large aniso- 
tropy we have extended their calculation, taking the first five K,, into account, 
and we made a comparison for some crystals with the result obtained when 
only the first three or four K, are considered. 

When five Ky are used v;(0¢) must be evaluated from (3) for five different 
directions, for four K,, four directions and for three K,, three directions. If n 
approaches one, % in (4) becomes independent of the choice of directions and — 


* The expression for K, contains a misprint. 
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of the number of terms in the expansion, but when yn is far from one, a will 
depend on the number of terms as well on the directions chosen in so far as 


the number of K,, used is not sufficient. 


On account of this % has been cal- 


culated for different combinations of directions, when three and four K,, are 


considered. For three K,, 
three combinations, which are in Table 1. 


five combinations have been used and for four K,, 


Table 1. The directions used in the calculation of 6p. 


Cn = Zn (C41 — Cqa) + Cay 
+2(n— 
Ly mpl Va b2 soout ° 
3 
C4 


p= —h(1+3K°) 
q= —ig (143K? +4 K°) 


1) | for n=1, 2,3 
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For the sake of brevity the expression for J is given only in the case of 
five K,,. 
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Table 2 contains the values of 0p obtained by this method. 


Table 2. The results for 6p. 


3 Km 4Km 5 Km 

” | I | II | III | IV | V VI | VII | vier 1x 

Na at 80° K 8.2 | 166 | 137 | 160 | 159 | 128 || 155 | 140 | 143 || 146 
Na at 290° K" | 7.0 | 245 | 123° | 141 }o24))| 116/138 (9126, ]220. saan 
K. . .« ».| 68c} 100} 86 | <98 | 98 } Si@pus6, | 88.) ou mmmeee 
Cu, y » | 32 | 341 | 326 | 334 | 338 | 316 || 336.2] 328.5] 330.6] 331.9 
Au 3; » » | 29 | 158° | 185:) 158 | 16% | 151; 160.0) 156.7 /°157-6))| Tose 
Ag . w» » | 29 | 217 | 209 | 213 | 216 | 204 || 214.6} 210.5] 211.5]) 212.3 
Al 4, os 9 | 22 | 407 | 407 | 407.,| 407) | 407 4} 407.2 | 407.2] 407.2) 2072 
NaCl ,, ,, , | 0.72 | 303 | 302 | 302 | 303 | 302 || 302.7] 302.7| 302.7|| 302.7 
NaCl at 0° K 0.56 | 322 | 321 | 318 | 321 | 319 || 321.2] 321.7] 321.6|| 321.5 
KCl at 280° K 0.36 | 224 | 224 | 215 | 224 | 219 || 223.3] 225.8] 225.1 | 224.6 
KCl at 0° K 0.29 | 230 | 230 | 216 | 230 | 222 || 228.9] 233.1] 232.0]| 231.2 


The elastic constants are taken from Hearmon [2] except those of NaCl, 
which are taken from Overtone and Swim [6]. They are all measured at room 
temperature except those of Na, KCl and NaCl. For Na the elastic constants 
at 80° K and 290° K have been used, for KCl and for NaCl those at 280° and 
290° K respectively have been used as well as those extrapolated to 0° K, for 
NaCl by Overtone and Swim and for KCl by Durand. In Table 2 one can see 
that combination II in Table 1 gives the best result, when only three K,, are 
used. When four K,, are used, the results are more stabilized, and thus the 
result for five K,, should be rather reliable. From Table 2 it is clear that the 
use of only three K,, is a good approximation for small anisotropy. 


The method of Hopf and Lechner and related methods 


The expression (2) can be written as follows: 


40 
% 1 dQ 
I=4 e ph Ne tee 
* eae J ~ (C+u)” 42 (6) 


Cc 
where rp pe Bs ies 
Cy — Cua 


From this Hopf and Lechner [3] computed the mean Debye characteristic 
temperature by approximating the integrand in (6) by means of a polynomial 
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of fifth degree. They evaluated the coefficients in the polynomial with La- 
ee of interpolation for O0<w<1. After that they only had to cal- 


culate i u; (n=1...5) from the coefficients in (3). This method has recently 


been tener by ne and Sutton [7]. These authors also used a poly- 
nomial of fifth degree but chosen so as to be exact or nearly exact at the 
points where the distribution functions of the roots have their maxima. In a 
later paper [9] Sutton has computed a correction term to (6) and has given this 
as a function of K and C on a arp He has considered the range where 
C>0.25 and where —0.85<(1—K)/2C<0 or 0<(1—K)/2<0.5. This range 
covers many cubic crystals as for example, Na, K and KCl. Outside the 
considered range the approximation with a fifth degree polynomial is not so 
good and, therefore, this method is then difficult to use. 

Schréder and v. Hagenow [8] have written Quimby and Sutton’s formula as 


w | Ze aa 
i \4aAJ | 


ss J(C, K) (7) 
0 


22k 
where Op = B gage 
h 


Op. 

They have evaluated J for different values of C and K in the region where 
0=C<2.8 and 0<K<4.5. The values obtained for J are written in a diagram. 
When C and K are known, J can easily be found from the diagram, and 0p 
can be computed from (7). The accuracy is, however, only 2 %. The considered 
range contains several crystals but for example Na lies outside, when Hearmon’s 
data for the elastic constants are used. 

In our investigation the correction term to Quimby and Sutton’s formula has 
been computed by expanding in K,, the difference between the integrand in (6) 
and the fifth degree polynomial. This has been done by means of three K,,, 
and 9p was calculated for Na and KCl. The results were 146° for Na (i.e. the 
same as with formula (5)) and 232.0° for KCl at 0° K. 

In order to compare the method of Quimby and Sutton with formula (5) 6p 
has been calculated for Au with exactly the same constants. The results were 
the same for the two methods; 0)=158.2°. 

For KCl at 0°K the value found for 6p is in ge agreement with that 
found by experimental measurements of C, below 2) K. Keesom and Pearl- 
man [4] found §,=233+3° K. They calculated 6) from the T°®-law, which 
must give a very good approximation in this region. 

Formula (5) can be used for every cubic crystal but the convergence is, of 
course, better the closer 7 is to unity. The error introduced by this method 
is, however, insignificant in- view of the uncertainty in the experimental deter- 
mination of the elastic constants, even when 7 is rather far from unity. 
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